Cuprate core level line shapes for different Cu-O networks 
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We have studied the Cu core level photoemission spectra in the Anderson impurity model for 
cuprates with different Cu-O networks, dimensionalities (zero, one, two and three) and Cu valencies 
(two and three). We focus on the shape of the leading peak and obtain very good agreement with 
the experimental data. We show how the shape of the spectrum is related to the valence electronic 
structure and the Cu-O network but also that other atoms can play a role. 



Boske et al. have shown that for a number of Cu 
compounds the shape of the core level photoemission 
spectrum strongly depends on the type of Cu-O network. 
This illustrates that it is not sufficient to think of the core 
level spectra as primarily atomic-like, but that there is a 
substantial influence from the valence electron structure. 
The core level spectrum can therefore give interesting in- 
formation about the electronic structure. 

The Cu compounds studied by Boske et al. are diva- 
lent, with Cu having mainly d 9 character. When a Cu 
core hole is created, it is energetically favorable to trans- 
fer an electron from the surrounding to the site of the core 
hole, forming a d 10 state. This leads to a spectrum with 
a leading peak of primarily d 10 character and a satellite 
of mainly d 9 character [Q . The main peak corresponds to 
a process where an electron has hopped from the valence 
band to the site with a core hole, and its shape reflects 
the valence band structure HQ . This leads to a broaden- 
ing of the main peak and possibly to further structures. 
Early experiments with a moderate resolution showed lit- 
tle structure. The high-resolution experiments of Boske 
et al. [Q, however, reveal a clear system dependence of 
the core level spectra. In the present paper we calculate 
these spectra and show that they are in good agreement 
with experiment. We explain in detail the origin of the 
differences between different cuprates. It has been sug- 
gested that the shape of the leading peak is explained in 
terms of a two-peak structure, due to so-called local and 
nonlocal screening H. The recent experiments, however, 
show that the shape can be substantially more compli- 
cated. 

To calculate the core spectrum, we use the Anderson 
impurity model 
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where v is a combined spin and orbital index. The (host) 
valence states are labelled by the energy e and v, and the 
Cu impurity 3d states are labelled by v p). The impu- 
rity level is pulled down from Sd to Ed — Udc when a core 



hole is created. The electrons on the impurity have the 
Coulomb repulsion U. 

It is crucial to use realistic hopping matrix elements 
V v (e). The local density of states (LDOS) p v (e) of a 3c? 
orbital with the label v is calculated for all compounds 
studied here, using the the local density approximation 
(LDA) and the LMTO method §. Then § 



J e — e — ^0 



(2) 



An Anderson impurity calculation using these matrix el- 
ements contains effects of all the 3c? orbitals, not just 
the ones on the impurity. An impurity calculation with 
U = exactly reproduces the LDOS of a lattice band 
structure calculation ||. In our calculation with U ^ 0, 
the many-body effects are treated explicitly on the site 
where the core hole is created, while it is assumed that 
the 3d states on the other sites can be treated in the LDA. 
In the discussion of the Anderson impurity model, it is 
often overlooked that all 3c? states are (approximately) 
included. 

We set U — 8 eV and Ud c — 10 eV and adjust Ed so 
that the relative weights of the main peak and the satel- 
lite agree with experiment. We have neglected multiplct 
effects. These effects are important for the shape of the 
satellite, but not for the shape of the main peak, which 
we consider here. 

To solve the model, we expand the ground-state as 



\E ) = A{^ v \d W ) 



dea(e)^|c/ lu )), 



(3) 



where the parameters A and a(e) are determined varia- 
tionally. |c? 10 ) is a state with all the valence states (|e^)) 
up to the Fermi energy (sjr) as well as all the impurity 
states filled. The first term in (||) corresponds to a d 9 
configuration and the second to a state where a valence 
electron with the energy e has filled the 3c?-hole, in the 
following referred to as d 10 e~ 1 . Then 



a v {e) = 



Vu(e) 



AE -e d + e' 



(4) 



AE is the lowering of the total energy due to hopping rel- 
ative to the energy of the d 9 state. Typically \AE — ea\ 
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is of the order of a few eV. A lowering of e leads to a de- 
crease of a(e) 2 , unless the variation of the denominator is 
compensated by a strong energy dependence of |y(e)| 2 . 
By adding more states to the Ansatz (||), the model can 
be solved increasingly accurately. However, the effect of 
additional states is essentially just a renormalization of 
Sd, which we treat as an adjustable parameter anyhow. 
The core level spectrum is calculated, by using a resol- 
vent operator and by inserting intermediate states simi- 
lar to the two terms in Eq. (||) . The spectra were 
broadened by a Lorentzian (FWHM=1.1 eV) to simulate 
life-time effects and by a Gaussian (FWHM=0.4 eV) to 
simulate the instrumental broadening. 

The results are shown in Fig. |l|, and they are in very 
good agreement with the experiments of Boske et al. 
0. In agreement with experiment, Ba3Cu204Cl2 and 
Li2Cu02 have very narrow peaks, Bi2CuC>4 a slightly 
broader peak and the remaining compounds significantly 
broader spectra. The shoulder of SrCu02 is more pro- 
nounced than for Sr2Cu03, as it should, and the three 
peak structure of Ba2Cu304Ci2 is reproduced very well. 
The main difference to experiment is that the leading 
structure of S^CuO^C^ is narrower. The high binding 
energy structures of S^CuOa and SrCuC>2 furthermore 
occur at somewhat lower binding energy and the high 
binding energy structure of SrCuC>2 is somewhat more 
pronounced. 

We now discuss the shape of the spectra. In the pres- 
ence of a core hole, the d 10 e~ 1 configuration is signifi- 
cantly lower in energy than the d 9 configuration. For the 
qualitative discussion (but not in the calculations) we as- 
sume that the coupling between these configurations can 
be neglected in the final state. Although this underesti- 
mates the weight of the leading peak, it reproduces its 
shape quite well (compare 'Spectrum' and l a 2 (e)' in Fig. 
||). The corresponding spectrum is 



p{e) = A 2 [a 2 (e)Q(-e) + S(e -e d + U dc )}, 



(5) 



where the first term gives the main peak and the second 
the satellite. The shape of the main peak is then deter- 
mined by a 2 (e), which is related to |l^(e)| 2 via Eq. (|4|). 
We therefore first study the relation between V(e) and 
the electronic structure. 

V(e) gives the hopping between the impurity and host 
(the rest of the solid but with the impurity removed). 
We first consider the case when the LDOS p»(s) in Eq. 
(0) is given by a semi-ellipse 
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where 2B is the width of the band. Performing the 
Hilbert transform in Eq. (|J) then gives the same shape 



ir\V v (e)\ 



Z7T 



B^ 



The overall strength of the hopping is determined by the 
width of the LDOS. This is expected, since a large width 
of the LDOS implies a strong hopping. 

The most important energy bands in the cuprates stud- 
ied here consist of antibonding combinations of one Cu 
orbital (e.g., x 2 — y 2 ) with p orbitals on the neighbor- 
ing O sites located around Ef and corresponding bond- 
ing combinations some 5-6 eV below Bp. The bonding- 
antibonding splitting is mainly determined by the Cu-0 
hopping, while the width of the two bands mainly de- 
pends on the hopping to more distant CuO„ units. As a 
model for the LDOS we use two semi-ellipses 
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FIG. 1. a) Theoretical core level spectra for different di- 
valent and one trivalent (LaCuOs) cuprate. The figure only 
shows the main peak, while the satellite at higher binding 
energies is not shown, b) Experimental spectra by Boske et 
al. Ilf for the divalent Cu compounds and by Mizokawa et al. 
[jn] for the trivalent LaCuOs. 
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\e\ < B. (7) 



This LDOS may result in a (5-function between the two 
continuous parts of the LDOS. If a|eo| >> B and 
(1 — a)|eo| » B, there is a (5-function at e = aso with 
the weight 7ro;(l — q)sq. This corresponds to the coupling 
to a linear combination of p orbitals on the neighboring O 
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atoms. The position of the 5-function depends, however, 
on the weights and widths of the bonding- and antibond- 
ing bands. In addition, there is a continuous part at the 
bonding and antibonding bands with a strength that is 
related to the widths of these bands. In practice, the 
LDOS is not zero between the bonding- and antibonding 
bands, and the ^-function then has a broadening. 

In Fig. we show results for Ba2Cu3 04Cl2, which has 
two inequivalent Cu atoms (Cu,4 and Cus). One of the 
Cu atoms (Cub) is connected to neighboring Cus atoms 
only via Cu-O-O-Cu bond and to neighboring Cu^ atoms 
via 90° Cu-O-Cu bond angels. Neighboring Cu atoms 
then couple to different O p orbitals, and they only have 
a weak coupling via other orbitals. The bonding- and an- 
tibonding bands are therefore narrow, as illustrated by 
the panel LDOS in Fig. |a. The resulting \V(e)\ 2 has a 
broadened (5-function between the bands but very little 
weight at the energies of the bands. The corresponding 
a 2 (e) also has most of the weight between the two bands. 
This shape agrees well with the result of the full calcu- 
lation (without the approximation (0)), and explains the 
narrow peak in the core level spectrum from Cug. The 
other Cu atom (Cua) is connected to the neighboring 
Cu atoms via 180° Cu-O-Cu bond angels. Neighboring 
Cu atoms then couple to the same O p orbital, i.e., they 
have a strong coupling. As a result, in particular the an- 
tibonding band is broad, as illustrated by the LDOS in 
Fig. ||b. The corresponding |y(e)| 2 has a broadened 6- 
function, which is shifted towards lower energies than for 
Cub, due to the different widths and weights of the two 
bands. In particular, however, there is a substantial cou- 
pling in the energy range of the antibonding band. This 
coupling is strongly enhanced by the energy denomina- 
tor (AE — Ed + e) of Eq. (^) in the calculation of a(e). 
Therefore, this coefficient is rather large for small bind- 
ing energies. Due to a structure in |T^(e)| 2 at s ~ — 2 
eV, a 2 (e) also has a structure at this energy. Finally, the 
large structure in |y(e)| 2 at e ~ —4 eV leads to a third 
structure in the spectrum. Adding the spectra from Cua 
and Cub gives the spectrum in Fig. ^. 

Following the analysis in Fig. we divide the sys- 
tems in roughly two classes. One class consist of systems 
where the CUO4 units have a weak coupling to the other 
CUO4 units, either due to the Cu-O-Cu bond angels being 
(rather close to) 90° or because the systems is approx- 
imately zero-dimensional (Od). In the second class the 
CUO4 units have a strong coupling to other CUO4 units, 
due to Cu-O-Cu bond angels of the order of 180°. 

Bi 2 Cu0 4 (Od), Ba 3 Cu20 4 Cl2 and Li 2 Cu0 2 belong to 
the first class. We would then expect a narrow peak in 
the core spectrum, occurring at a higher binding energy 
than the leading edge. This is confirmed for the latter 
two systems, while both theory and experiment shows 
a somewhat broader peak for Bi 2 Cu04. The reason for 
this is that there is a coupling between the CUO4 units 
via the Bi atoms. As a result the structure in |F(e)| 2 



between the bonding- and antibonding bands is broader, 
leading to a somewhat broader core level spectrum. This 
illustrates that not only the Cu-0 network plays a role. 

Sr 2 Cu0 2 Cl 2 , Sr 2 Cu0 3 and SrCu0 2 belong to the 
second class. In the two-dimensional (2d) structures 
Ba 2 Cu304Cl 2 (Cua) and Sr 2 Cu0 2 Cl 2 the antibonding 
band is broader than in the Id structures Sr 2 Cu03 and 
SrCu0 2 , due to the more efficient hopping to neighboring 
CUO4 units in the 2d case. Therefore, the structure in 
V^(e)| 2 at e ~ — 2 eV is broader and the peak in |F(e)| 2 
at e ~ —4 eV has been pushed to somewhat lower ener- 
gies in the 2d case. The large binding energy structure in 
the main peak is then pushed towards somewhat larger 
energies in the 2d case, in agreement with experiment 
(see C and C' in Fig. El). 
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FIG. 2. Core level spectrum, broadened a 2 (e), [F(e)| 2 and 
local density of states (LDOS) for the two ( a) and b)) in- 
equivalent Cu atoms in Ba2Cu3C>4Cl2. |V(e)| 2 and the LDOS 
have been given a small (FWHM=0.03 eV) broadening, while 
a 2 (e) has been given the same broadening as the spectrum. 
We have put sf — 0. 



Above we have considered divalent Cu. We have earlier 
studied the trivalent NaCu0 2 , and found that this has a 
narrow leading peak due to a narrow structure in |I^(e)| 2 
pj. An alternative explanation was suggested by a clus- 
ter calculation, where it was argued that the reason for 
the small width is the valence (3+) ||. We have there- 
fore considered a different trivalent compound, namely 
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LaCu03- The corresponding spectrum is shown at the 
bottom of Fig. Q. The main peak is broad, both accord- 
ing to theory and experiment J5~T[ ] . The reason is similar 
as for the broad peaks in the divalent compounds, namely 
a broad antibonding band due to a Cu-O-Cu bonding an- 
gle of about 180°. Due to the three-dimensional charac- 
ter of this compound, the band is unusually broad, but 
only a relatively small fraction is occupied, because the 
compound is trivalent. Therefore the broadening of the 
leading peak is only modest. If we had instead studied a 
model of the type used in the cluster calculations, with 
only one orbital per atom, the antibonding band would 
have been empty. The Anderson impurity model would 
then have predicted a narrow leading peak, as the clus- 
ter calculations do. To obtain the correct description 
of LaCuC>3 it is therefore crucial to use a realistic model 
which includes the partial filling of the antibonding band, 
and which goes beyond the normally used cluster models. 

In an exact diagonalization treatment of a CU3O10 
cluster, the leading peak consists of two ((-functions. 
These correspond to so-called local and nonlocal screen- 
ing II , indicating that the charge screening of the core 
hole comes from the closest O4 unit or from the neighbor- 
ing CUO4 units, respectively. Similar results are obtained 
in the Anderson impurity model for such a cluster. Find- 
ing the solutions of the host, we notice that the highest 
occupied level has a strong antibonding character and is 
located on the neighboring CUO4 units. The second high- 
est occupied level (with a coupling to the impurity) has 
a large weight on the O4 unit around the impurity. Due 
to the missing hybridization with the impurity, this host 
level is less antibonding and lower in energy. The two 
^-functions in the spectrum essentially correspond to an 
electron from one of these two host levels hopping into 
the impurity orbital. Due to their different characters, 
these levels give local or nonlocal screening. 

It has been emphasized that it is important to take 
into account the formation of a Zhang-Rice singlet on the 
CUO4 units neighboring the unit with a core hole @ . We 
find that in the Anderson impurity model of a CU3O10 
cluster, the two lowest final states contain comparable 
amounts of a Zhang-Rice singlet as is found in exact di- 
agonalization for such a cluster. If the highest occupied 
host level is emptied, a large amount of antibonding char- 
acter is removed from the neighboring CUO4 units. This 
alone is enough to give a large overlap to a state with a 
Zhang-Rice singlet on one of these units. If instead the 
second highest host level is emptied, the screening charge 
is primarily taken from the local O4 unit, and the local 
singlet character on the neighboring CUO4 units remains 
small. 

To further test the role of the singlet in core level pho- 
toemission, we have performed a calculation for a CU3O10 
cluster with 25 valence electrons. This system has only 
one valence hole, and a Zhang-Rice singlet cannot form. 
The change of the number of electrons gives a somewhat 



different spectrum, but the main peak still has a substan- 
tial broadening, illustrating that just going beyond a sin- 
gle site model (CUO4) leads to a broadening of the leading 
peak. The considerations above illustrate the usefulness 
of a picture where an electron from a host level hopps 
into the impurity orbital. 

It has been found that even rather 'large' clusters like 
CU3O10 may not be converged with respect to their size 
Jjjfl . The inclusion of additional orbitals like O may 
also be important fL2f . Here we have found further exam- 
ples of a need for more realistic models, e.g., for Bi2CuC>4 
and LaCuC>3. In an exact diagonalization approach, how- 
ever, it is not possible to make the cluster much larger 
or to add many other orbitals, since the many-electron 
space then becomes too large, while such a system can 
be treated in the Anderson impurity approach. 

To summarize, we have calculated core level spectra for 
a number of cuprates in the Anderson impurity model 
using ab initio hopping matrix elements. In this way, 
many-body effects are treated explicitly for the Cu site 
where the core hole is created, while the other Cu sites 
are treated in the LDA. This model gives a very good 
agreement with experimental core level spectra. The 
shape of the leading peak is influenced substantially by 
the valence electronic structure and thereby the lattice 
structure. The presence of (approximately) 180° Cu-O- 
Cu bond angels leads to a broad and structured lead- 
ing peak, and an increase in the dimension of the Cu-0 
network tends to increase the width. In particular for 
the trivalent compounds, it is important to go beyond 
the simple cluster models with one orbital per site, since 
these models always give a narrow leading peak. 

This work has been supported by the Max-Planck- 
Forschungspreis . 
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